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John  A.  Guinn,  Jr.,  Gilbert  N.  Plass,  and  George  W.  Kattawar 


Abstract 

Time-averaged  intensities  are  computed  for  the  glitter  pattern  of  sun- 
light on  a wind-ruffled  sea.  In  the  present  work  the  azimuthal ly  asymmetric 
form  of  the  wave-slope  distribution  measured  by  Cox  and  Munk  is  used,  so  that 
the  glitter  pattern  depends  on  the  wind  direction  as  well  as  magnitude. 

Results  are  calculated  for  both  a clean  ocean  surface  and  one  covered  with 

an  oil  film.  The  glitter  pattern  is  quite  different  with  the  oil  film  due 
to  the  quieting  action  of  the  oil  on  the  waves,  the  greater  reflectivity  of 
the  surface  when  oil  is  on  the  surface,  and  the  larger  values  of  polarization 

for  the  reflected  light  at  most  typical  angles  in  the  glitter  pattern.  These 

effects  are  illustrated  and  discussed  as  possible  methods  for  the  detection 
of  oil  slicks.  The  glitter  pattern  observed  at  satellite  altitudes  is  also 
calculated. 
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I.  Introduction 

Glitter,  the  pattern  of  coruscations  one  can  see  on  a wind-ruffled  water 
surface  that  is  illuminated  by  a strong  source,  is  a commonplace  phenonmenon 
of  the  type  that  M.  Minnaert  particularly  enjoyed  investigating.  Although 
frequently  observed  and  often  recorded  --  sometimes  inaccurately  --  by  artists, 
little  thought  has  been  devoted  to  glitter  pattern  morphology  and  its  causes. 
The  existence  of  a glitter  pattern  on  a sunlit  sea  is  indeed  not  surprising, 
but  it  is  not  entirely  obvious  why  the  pattern  should  assume  a particular 
configuration  under  a given  set  of  circumstances.  Through  some  basic  geo- 
metrical considerations,  Hulburt^  related  the  solar  elevation  angle  and  what 

he  called  the  effective  maximum  wave  slope,  a measure  of  sea  roughness,  to 
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the  width  of  the  glitter  pattern.  Minnaert  ’ added  considerable  detail  to 

the  analysis  of  the  problem;  his  observations  are  discussed  further  as  they 

apply  to  these  results.  Certain  geometrical  insights  regarding  glitter  and 

related  comments  pertaining  to  the  distortion  of  the  shapes  of  objects  re- 

fleeted  in  a wind-ruffled  water  surface  have  been  provided  by  Tricker  and 
5 6 

Goodell  . Gambling  has  studied  sun  glitter  in  the  infrared  spectral  region. 

Plass  et  have  calculated  the  time-averaged  intensities  for  the 
glitter  pattern  of  sunlight  on  a wind-ruffled  sea.  Isopleths  drawn  on  graphs 
simulate  glitter-pattern  photographs.  The  pattern  is  computed  from  the  wave- 
surface  orientation  required  for  direct  reflection  from  the  source  to  the 
observer;  the  intensity  is  computed  from  the  probability  of  occurrence  of 
this  orientation,  as  determined  from  the  Cox-Munk  distribution,  together  with 
the  fraction  of  the  radiation  reflected  at  that  particular  angle  of  incidence. 
The  curvature  of  the  earth  is  taken  into  account. 
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II.  Method  of  Calculation 

The  glitter  is  composed  of  the  many  bright,  instantaneous  points  of 
light  created  on  the  water  surface  by  local  surface  orientations  which  are 
momentarily  appropriate  for  the  direct  reflection  of  the  source  into  the  eye 
of  the  observer.  The  brightness  of  each  glint  depends  on  the  reflection 
coefficient  of  the  reflection  which  produces  the  glint;  the  reflection  co- 
efficient in  turn  depends  on  the  angle  of  incidence  of  the  reflection.  The 
perceived  (time-averaged)  brightness  of  a particular  region  of  the  glitter 
pattern  is  a consequence  of  the  brightness  of  individual  glints  and  their 
frequency  of  occurrence.  The  frequency  of  occurrence  of  highlights  at  a 
particular  location  on  the  water  surface  depends  on  the  local  surface 
orientation  required  to  produce  the  highlights  and  the  nature  of  the  pre- 
vailing distribution  of  surface  orientations.  A glitter-pattern  model  there- 
fore must  assume  a surface-orientation  distribution. 

In  the  present  work  the  general  procedure  described  in  Ref.  7 is  followed, 
but  the  new  calculations  take  into  account  the  complete,  azimuthally  asymmetric 

g 

form  of  the  wave-slope  distribution  measured  by  Cox  and  Munk  . A bivariate 
normal  distribution  having  different  standard  deviations  for  the  upwind  (o^) 
and  crosswind  (a^)  components  is  used.  This  distribution  is  modified,  to 
produce  a Gram-Charlier  distribution,  by  multiplying  the  bivariate  normal 
density  function  by  a sum  of  products  of  Hermite  polynomials.  The  density 
function  is  therefore 

p{2c.Zu)  = [G(a,b)/2TTo^a^^]  exp[-%(a^  + b^)].  (1) 

In  Eq.  (1),  z and  z are  the  crosswind  and  upwind  components,  respectively, 

w U 

of  the  wave  slope,  a = and  b = z^/Oy*  6(a,b),  a sum  of  products  of 
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Hermite  polynomials,  Is  the  Gram-Charlier  sum,  which  can  be  written 

G(a,b)  = 1 -%C2i(d^-  l)b-CQ3(b^-  3b)/6  + c^gCa^  - 6a^  + 3)/24 

+ %C22{a^-l)(b^-l)  + CQ4(b‘^-6b^  + 3)/24  + . . . (2) 

Notice  that  two  of  the  terms  are  asymmetric  about  the  origin  in  the  upwind 
direction.  The  coefficients  of  the  sum,  experimentally  determined  by  Cox 

O 

and  Munk  , are  given  in  Table  I for  a clean  water  surface  and  for  a surface 
covered  by  an  oil  slick. 

The  absence  of  a cross-product  term  in  the  argument  of  the  density- 
function  exponential  is  due  to  the  choice  of  the  upwind  and  crosswind  direc- 
tions as  the  coordinate  axes.  There  is,  in  other  words,  no  statistical 
correlation,  in  the  absence  of  the  Gram-Charlier  sum,  between  the  upwind 
and  crosswind  components  of  the  wave  orientations. 

Q Q 

The  data  of  Cox  and  Munk  (as  well  as  those  of  Duntley  who  achieved 
similar  results)  indicate  a linear  relationship  between  the  upwind  and  cross- 
wind  variances  and  the  wind  speed.  The  Cox  and  Munk  results  for  a clean 
water  surface  are 

o ^ = 0.003  + 0.00192W 
c 

and  (3) 

= 0.000  + 0.00316W, 

while  for  a surface  covered  by  an  oil  slick 

= 0.003  + 0.00084W 

and  (^) 

2 


= 0.005  + 0.00078W, 


where  W is  measured  in  m/sec  at  a height  of  12,5  m above  the  sea  surface. 

For  the  clean  surface  the  crosswind  constant  term  dominates  for 
W < 2.42  m/sec,  while  ^ fof  greater  wind  speeds.  For  the  slick  surface 
Ou > Oc  for  W < 33.3  m/sec,  which  completely  covers  the  range  of  wind  speeds 
considered  here.  Thus  the  asymmetries  have  the  same  sense  for  2.42  m/sec  < 

W < 33.3  m/sec  and  the  opposite  sense  for  W < 2.42  m/sec.  However,  Eqns.  (3) 
and  (4)  show  that  the  upwind-crosswind  asymmetry  is  very  small  for  the  slick 
surface  compared  to  the  clean  surface.  The  equations  also  show  the  smoothing 
effect  of  the  oil  slick. 

The  graphs  presented  here  represent  the  projection  of  the  sea  surface 
onto  an  inclined  plane  placed  in  front  of  the  observer  and  simulate  the  image 
that  would  be  recorded  on  a photographic  plate.  Each  point  on  the  plane 
corresponds  uniquely  to  a point  on  the  sea  surface  (except  for  points  which 
project  above  the  horizon).  The  correspondence  is  determined  by  requiring 
the  lines  of  projection  to  meet  at  the  observer.  Calculations  were  made  of 
glitter-pattern  intensity  (time-averaged),  polarization,  and  certain  other 
quantities  for  points  in  a rectangular  grid  on  the  inclined  plane.  For  each 
of  these  points  the  corresponding  sea-surface  location  was  determined,  and 
from  that  the  surface  orientation  required  for  a direct  reflection  of  the 
source  to  the  observer  calculated.  The  required  surface  inclination  and 
azimuth  were  then  applied  to  the  density  function  already  discussed  to 
calculate  the  probability  of  occurrence  of  the  surface  orientation  in 
question.  This  probability  then  became  a factor  in  the  calculation  of  the 
glitter  intensity  at  the  grid  point,  as  did  the  coefficient  of  reflection 
determined  from  the  angle  of  incidence  of  the  light  from  the  source  on  the 
surfiico  rlemont  of  correct  oriental  ion.  The  computed  quantities  assigned 
to  the  grid  points  were  used  to  geni'rate  the  contour  plots  shown. 


-6- 


The  coordinate  system  used  for  the  graphs  is  the  same  as  that  used  in 
Ref.  7,  which  provided  a detailed  description.  Very  briefly,  the  graphs  re- 
present a projection  of  the  glitter  pattern  and  a sea-surface  rectangular 
coordinate  system  onto  a plane  inclined  so  that  it  is  perpendicular  to  the 
line  connecting  the  observer  with  the  center  of  the  vertical  axis  of  the 
graph.  Since  these  graphs  generally  are  not  bilaterally  symmetric,  as  were 
those  of  Ref.  7,  both  sides  of  the  principal  plane  are  represented.  Wind 
direction  is  measured  clockwise  from  the  solar  azimuth. 

For  a clean  water  surface  the  coefficient  of  reflection,  for  a single 
interface,  is  simply  determined  from  the  angle  of  incidence  and  the  index 
of  refraction  of  sea  water  (assumed  to  be  1.338).  The  coefficient  of  re- 
flection is  slightly  more  complicated  when  there  is  a double  interface,  as 
is  the  case  when  the  sea  surface  is  contaminated  with  an  oil  slick.  Let 
R^  and  R2  represent  the  fractional  intensities  reflected  at  the  first  and 
second  interfaces,  respectively.  If  the  interfaces  are  assumed  parallel,  a 
simple  calculation  shows  that  the  total  reflected  intensity  including  all 
orders  of  multiple  reflection  within  the  oil  film  is 

R = R^  + R2(1  -R^)^  / (1  -R1R2).  (5) 

Interference  effects  and  absorption  and  scattering  within  the  layer  are 
neglected. 

The  coefficient  of  reflection  for  a single  interface  (air-water)  is 
shown  in  Fig.  1 (solid  curve);  for  a double  interface  (air-layer-water,  where 
the  layer,  representing  an  oil  slick,  has  an  index  of  refraction  of  1.5)  it 
is  shown  as  a broken  line  with  solid  circles.  A similar  double  layer  with 
oil  with  index  of  refraction  of  1.6  is  shown  as  a dashed  line.  The  parallel 
component  for  water  alone  is  zero  at  the  Brewster  angle,  53.23°,  but  with  a 
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double  interface  the  parallel  component  never  quite  reaches  zero,  as  shown 
in  the  inset  on  the  graph.  For  oil  with  n = 1.5,  the  minimum  value  for  the 
parallel  component  occurs  at  an  angle  of  incidence  of  56.6°  and  the  maximum 
polarization  is  0.9909.  When  n = 1.6,  the  corresponding  angle  is  58.5°  and 
the  maximum  polarization  is  0.9824.  The  curve  for  the  latter  case  is  not 
shown  in  the  inset  as  it  is  off  scale. 

The  polarization  of  the  reflected  radiation  for  the  same  three  cases 
is  shown  in  Fig.  2 with  the  assumption  that  the  source  is  unpolarized. 

III.  Intensities  for  Observer  Near  Surface 

The  variation  in  the  structure  of  the  glitter  pattern  with  changes  in 

the  wind  direction  are  illustrated  in  this  section.  Figs.  3,  4,  and  5 are 

parametrically  identical  except  for  wind  direction.  Each  figure  represents 

a clean  surface  illuminated  by  a source  (at  infinite  distance,  i.e.  the  sun) 

at  30°  elevation,  ruffled  by  a 10  m/sec  wind,  and  observed  from  a height  of 

5 m above  the  sea  surface.  The  wind  directions  are  0°,  45°,  and  90°, 

respectively.  The  principal  plane  corresponds  to  the  centerlines  of  these 

figures  and  wind  directions  are  measured  from  this  line.  The  coordinates 

0 and  e are  explained  in  Ref.  7.  The  numbers  associated  with  each  contour 
X y 

line  give  the  logarithm  to  the  base  10  of  the  intensity  of  the  time-averaged 
radiance. 

Though  the  range  of  intensity  is  essentially  the  same  in  each  case,  the 
effect  of  varying  wind  direction  can  be  clearly  seen  in  the  morphological 
differences  among  these  figures.  The  U-shaped  contours  of  Fig.  3 (wind  from 
the  solar  azimuth)  give  way  to  V-shaped  contours  in  Fig.  5 (direct  crosswind). 
The  brightest  contours  of  Fig.  5 are  much  closer  to  the  triangular  forms  noted 

3 

by  Minnaert  than  any  of  the  contours  generated  by  the  azimuthal ly  symmetric 
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wave-slope  distributions  of  Ref.  7.  Fig.  5 is  bilaterally  symmetric  because 
of  the  small  magnitude  of  Cq^,  the  coefficient  of  the  anti  symmetric  upwind 
terms  of  the  Gram-Charlier  sum,  but  the  contours  of  Fig.  4 are  skewed 
noticeably  to  the  right  by  a diagonal  wind.  Because  all  crosswind  terms 
are  symmetric.  Fig.  3 displays  the  bilateral  symmetry  characteristic  of  all 
cases  for  which  the  wind  is  from  the  solar  azimuth. 

The  slick-surface  counterparts  to  Figs.  3 and  5 are  Figs.  6 and  7, 
respectively.  All  parameters  correspond  to  the  earlier  figures,  except  that 
the  surface  is  assumed  to  be  covered  with  an  oil  layer  with  an  index  of  re- 
fraction of  1.5.  Waves  of  lesser  inclination  are  favored  when  the  slick  is 
present,  hence  the  contours  are  drawn  in,  and  the  brightest  spot  is  no 
longer  on  the  horizon.  The  intensity  at  the  horizon  is  not  greatly  diminished 
from  the  maximum,  particularly  in  the  case  of  Fig.  6;  the  intensity  difference 
is  only  21%  between  the  maximum  and  the  brightest  contour  that  meets  the 
horizon.  The  glitter  pattern  would  be  a narrow,  bright  streak  extending 
downward  from  the  horizon.  The  contrast  between  the  maximum  and  the  horizon 
is  somewhat  greater  in  Fig,  7,  the  intensity  difference  with  respect  to  the 
maximum  being  44%  at  the  brightest  contour  to  meet  the  horizon;  the  glitter 
pattern  is  noticeably  broader  than  the  one  in  Fig.  6. 

The  next  two  cases  were  calculated  for  a 30°  source  elevation,  an  observer 
at  a height  of  5 m,  and  a 1 m/sec  wind  from  a 45°  azimuth.  The  only  difference 
between  the  two  is  that  Fig.  8 represents  a clean  surface  and  Fig.  9 a slick 
surface.  The  opposite  skewing  of  the  isopleths,  despite  the  common  wind 
direction  for  the  two  graphs,  is  due  to  the  fact  that  at  1 m/sec  the  crosswind 
variance  dominates  for  the  clean  surface,  while  the  opposite  is  the  case  for 
the  slick  surface. 

An  oblique  wind  does  not  assure  bilateral  asymmetry.  Fig.  10  (5  m 
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observer  height,  5°  sun  elevation,  oil  slick  surface,  10  m/sec  wind  speed, 
and  45°  wind  azimuth)  is  very  nearly  synunetric  because  the  upwind  and  cross- 
wind  variances  are  approximately  the  same  for  this  wind  speed. 

IV.  High-Altitude  Observer 

The  next  seven  figures  represent  the  calculated  appearance  of  the  glitter 
pattern  to  an  observer  at  a height  of  150  km.  In  each  case  the  sun  is  30° 
above  the  horizontal  reference  plane.  The  horizon  can  be  seen  as  the  upper- 
most curve  in  each  figure.  A wind  speed  of  10  m/sec  was  used  for  Figs.  11-13. 
Figs.  11  and  12  represent  a clean  surface,  while  Fig.  13  is  for  an  oil  slick 
surface.  In  Fig.  11,  for  which  the  wind  direction  is  45°,  the  asymmetry  is 
quite  clear,  with  the  bright  spot  on  the  horizon  distinctly  to  the  right  of 
the  principal  plane.  A wind  shift  to  90°  (Fig.  12)  eliminates  the  asymmetry 
and,  because  the  crosswind  variance  is  less  than  the  upwind  variance,  moves 
the  bright  spot  away  from  the  horizon.  It  also  produces  a shield,  one  of 
the  distinctive  shapes  generated  by  certain  combinations  of  parameters.  In 
Fig.  13  the  wind  azimuth  is  still  90°,  but  an  oil  slick  has  been  added.  This 
change  shifts  the  bright  region  farther  below  the  horizon  and  splits  the 
brightest  spot  into  two  parts. 

Figs.  14-17  were  generated  assuming  a 1 m/sec  wind  speed,  a 150  km 
observer  height,  and  a 30°  solar  elevation.  The  oil  slick  model  was  used 
for  Figs.  16  and  17.  The  wind  azimuth  is  0°  for  Figs.  14  and  16  and  90°  for 
Figs.  15  and  17.  At  1 m/sec  the  crosswind  standard  deviation  dominates  for 
the  clean  surface,  while  for  the  surface  with  an  oil  slick  the  upwind  standard 
deviation  is  larger.  (The  ratio  of  the  larger  figure  to  the  smaller  is 
roughly  1.5  in  each  case.)  As  a consequence  the  clean-surface  contours  are 
broadest  when  the  wind  azimuth  is  0°,  while  contours  for  the  slick  surface 
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are  broadest  for  a 90°  wind  azimuth.  Thus,  the  heart-shape  of  Fig.  14 
(clean  surface)  becomes  the  tear  drop  of  Fig.  16  (oil  slick),  while  the 
tear  drop  of  Fig.  15  becomes  the  heart-shape  of  Fig.  17  (oil  slick). 

Figs.  18  and  19  (for  which  the  Gram-Charl ier  sum  was  not  used)  were 
calculated  for  a nearby  source.  In  both  cases  the  wind  speed  is  10  m/sec, 
the  wind  azimuth  is  45°,  the  observer  height  is  5 m,  the  source  is  horizontally 
displaced  100  m from  the  obsetver,  and  there  is  no  oil  slick.  The  source 
for  Fig.  18  is  at  a height  of  62.74  m,  and  that  for  Fig.  19  is  13.75  m high. 

The  higher  source  has  an  apparent  elevation  angle  of  30°  as  seen  from  the 
observer's  position,  and  the  apparent  elevation  angle  of  the  lower  source 
is  5°.  The  slanting  pattern  of  Fig.  18  is  somewhat  like  that  of  Fig.  18  of 
Minneart's  book  (Ref.  3,  p.  22),  which  depicts  the  glitter  pattern  of  a lamp 
across  a canal.  It  is  particularly  apparent  in  Fig.  19  that  the  intensity 
drops  off  very  rapidly  near  the  horizon.  This  is  due  to  very  high  wave 
inclinations  required  for  reflection  toward  the  observer  for  waves  on  the 
far  side  of  the  source. 

V.  Polarization 

The  source  is  assumed  to  be  unpolarized.  The  degree  of  polarization 
induced  by  reflection  is  a function  of  the  angle  of  incidence  at  each  gi  id 
point.  The  wind  speed  is  irrelevant  for  polarization  since  the  effect  depends 
only  on  the  angle  of  incidence  at  the  wave  facet  for  each  reflection.  The 
probability  of  occurrence  of  the  required  wave  slope  does  not  enter  into 
the  calculation  of  the  polarization,  although  the  radiation  obviously  could 
not  be  observed  if  this  probability  is  too  small.  Fig.  20  shows  the  isopleths 
of  degree  of  polarization  expressed  in  percent  for  a source  displaced  100  m 
horizontally  from  the  observer,  a source  height  of  62.74  m,  an  observer 
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height  of  5 m,  and  a clean  water  surface.  There  is  a large  region  in  which 
the  reflected  light  is  strongly  polarized  around  the  line  for  which  the 
observed  light  is  lOOi  polarized.  This  completely  polarized  light  is  re- 
flected from  the  water  surface  at  the  Brewster  angle  of  53.23°  (for  n = 1.338). 
When  the  sun  is  the  source  (at  an  effective  infinite  distance)  the  polarization 
drops  from  the  Brewster  angle  to  the  horizon  (see  Figs.  22  and  23  of  Ref.  7). 

On  the  other  hand,  when  the  source  is  at  a finite  distance,  the  polarization 
first  decreases  toward  the  horizon  and  then  increases  to  another  maximum, 
as  the  light  is  once  again  reflected  from  the  wave  facet  at  the  Brewster  angle. 
For  the  case  shown  in  Fig.  20,  the  first  Brewster  angle  is  observed  at  50.19° 
from  the  vertical;  this  light  is  reflected  from  wave  facets  making  an  angle 
of  3.04“ to  the  horizontal.  For  the  second  Brewster  angle  the  corresponding 
angles  are  86.28°  and  33.05°,  respectively. 

Figure  21  is  parametrically  the  same  as  Fig.  20  except  that  an  oil  slick 
has  been  added.  Several  changes  are  evident.  As  discussed  in  Sect.  II,  the 
maximum  polarization  for  an  oil  slick  with  n = 1.5  is  99.09%  for  an  angle  of 
incidence  of  56.6°.  A comparison  of  Figs.  20  and  21  shows  the  decreased  value 
of  the  polarization  at  the  maximum  when  an  oil  slick  is  present  as  well  as 
the  displacement  of  this  maximum  toward  the  horizon.  For  the  oil  slick  case, 
the  observer  must  look  at  57.44°  from  the  vertical  to  observe  the  light  re- 
flected at  the  first  Brewster  angle;  this  light  is  reflected  from  wave  facets 
making  an  angle  of  0.84°  to  the  horizontal.  For  the  second  Brewster  angle 
the  corresponding  angles  are  85.76°  and  29.15°,  respectively.  Another 
difference  between  the  two  figures  is  the  larger  value  of  the  polarization 
for  the  oil  slick  case  in  the  large  region  between  the  two  Brewster  angles. 

In  this  region  the  angle  of  incidence  of  the  light  on  the  wave  facets  is 
greater  than  the  Brewster  angle.  Fig.  2 shows  that  the  polarization  for  a 
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given  angle  of  incidence  is  greater  for  the  oil  slick  case  than  for  clean  i 

1 

water  for  this  range  of  angles.  I 

i 

VI.  Possible  Detection  of  Oil  Slicks 

Several  possible  methods  for  the  detection  of  oil  slicks  are  suggested 
by  these  results.  A measurement  of  the  polarization  of  the  reflected 
radiation  in  the  brightest  part  of  the  glitter  pattern  would  show  a slight 
decrease  in  the  maximum  polarization  at  the  Brewster  angle  (99.09%  for  oil 
slick  instead  of  100%  for  clean  water  in  the  example  considered  in  Sect.  V); 
a shift  of  the  Brewster  angle  toward  the  horizon  when  the  oil  slick  is  present 
(by  7.25°  in  the  example);  higher  polarization  values  in  the  broad  region 
from  near  the  Brewster  angle  to  angles  near  the  horizon  (polarization  about 
5%  higher  in  this  region  for  the  oil  slick  in  the  example).  The  second 
Brewster  angle  discussed  in  Sect.  V is  only  observed  for  sources  at  finite 
distances;  the  differences  between  a clean  and  oil  slick  surface  are  less 
in  this  case  and  would  be  more  difficult  to  observe.  The  contrast  in  the 
polarization  measurements  for  adjacent  regions  of  the  ocean  with  and  without 
an  oil  slick  should  be  relatively  easy  to  measure.  The  great  advantage  of 
measuring  the  polarization  is  that  the  value  observed  at  a given  angle  is 
independent  of  the  sea  state;  of  course  the  probability  of  a given  wave 
slope  must  be  large  enough  to  reflect  measurable  light  to  the  observer. 

Oil  slicks  can  also  be  observed  by  the  changes  in  the  intensity  of  the 
reflected  light  as  discussed  in  Sect.  Ill  and  IV.  The  main  effects  come  from 
an  increase  in  the  reflectivity  due  to  the  higher  index  of  refraction  of 
oil  than  of  water  and  to  the  decreased  probability  of  higher  wave  slopes  when 
an  oil  slick  is  present.  The  resulting  changes  in  the  intensity  of  the  re- 
flected light  in  the  glitter  pattern  as  shown  in  the  figures  suggest  many 
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possibilitles  for  observing  the  slicks. 

All  of  the  calculations  discussed  so  far  assume  that  geometrical  optics 
applies  to  this  problem.  Effects  due  to  the  wave  nature  of  light  also  can  be 
prominent  with  an  oil  film.  The  many  colors  in  the  light  reflected  form  such 
a film  are  well  known.  These  occur  as  the  light  at  a particular  wavelength 
reflected  from  the  top  surface  of  the  film  interferes  either  constructively 
or  destructively  with  the  light  reflected  from  the  oil-water  surface.  Light 
form  a tuneable  laser  reflected  from  an  oil  slick  should  show  maxima  and  minima 
as  the  wavelength  is  varied,  corresponding  to  these  interference  effects. 

Such  maxima  and  minima  would  be  a clear  signature  of  an  oil  slick,  as  they 
do  not  exist  for  a clean  water  surface. 

VII.  Conclusion 

The  results  of  the  glitter  calculations  show  how  the  pattern  varies  with 
wind  direction  as  well  as  wind  velocity.  Both  the  shape  of  the  glitter  pattern 
and  the  intensity  of  the  reflected  light  change  appreciably  when  an  oil  slick 
is  on  the  water  and  this  fact  can  be  used  to  locate  such  slicks.  The  polar- 
ization of  the  light  in  the  glitter  pattern  is  also  appreciably  changed  when 
an  oil  slick  is  on  the  surface.  Results  are  shown  for  an  observer  just  above 
the  ocean  surface  as  well  as  at  satellite  altitudes.  The  radiation  in  the 
glitter  pattern  is  an  important  component  of  the  upwelling  radiant  energy 
above  an  ocean. 

This  research  was  supported  by  the  Office  of  Naval  Research  through 
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Captions  for  Figures 

Reflectivity  as  a function  of  the  angle  of  incidence  (0)  for  a clean 
water  surface  (solid  curve);  oil  film  (n  = 1.5)  on  water  (dashed 
curve  with  solid  circles);  oil  film  (n  = 1.6)  on  water  (dashed 
curve).  Note  the  different  scales  for  the  left  and  right  half 
diagrams.  In  each  case  the  lower  curve  is  for  the  parallel  component, 
while  the  upper  curve  is  for  the  perpendicular  component.  The 
inset  shows  the  reflectivity  of  the  parallel  component  near  the 
Brewster  angle. 

The  polarization  of  the  reflected  radiation  as  a function  of  the 
angle  of  incidence  (0)  for  a clean  water  surface  (solid  curve);  oil 
film  (n  = 1.5)  on  water  (dashed  curve  with  solid  circles);  oil  film 
(n  = 1.6)  on  water  (dashed  curve). 

Intensity  isopleths  for  a source  elevation  of  30°,  an  observer  height 
of  5 m,  and  a wind  from  a 0°  azimuth  at  10  m/sec.  The  numbers  beside 
the  curves  give  the  logarithm  of  the  intensity  to  the  base  10. 

Maximum  log^Q  intensity  is  -4.018.  The  coordinates  are  explained 
in  the  text  and  in  Ref.  7. 

Intensity  isopleths  for  a source  elevation  of  30°,  an  observer  height 
of  5 m,  and  a wind  from  a 45°  azimuth  at  10  m/sec.  Maximum  log.|Q 
intensity  is  -4.102. 

Intensity  isopleths  for  a source  elevation  of  30°,  an  observer  height 
of  5 m,  and  a wind  from  a 90°  azimuth  at  10  m/sec.  Maximum  log^Q 
intensity  is  -4.215. 

The  same  as  Fig.  3,  except  that  an  oil  slick  is  added.  Maximum 
log^Q  intensity  is  -4.211. 
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7.  The  same  as  Fig.  5,  except  that  an  oil  slick  is  added.  Maximum 
log^Q  intensity  is  -4.250. 

8.  Intensity  isopleths  for  a clean  surface,  a source  elevation  of  30°, 
an  observer  height  of  5 m,  and  a wind  from  a 45°  azimuth  at  1 m/sec. 
Maximum  log^Q  intensity  is  -4.009. 

9.  The  same  as  Fig.  8,  except  that  an  oil  slick  is  added.  Maximum 
log-jQ  intensity  is  -3.911. 

10.  Intensity  isopleths  for  an  oil  slick  surface,  a source  elevation 
of  5°,  an  observer  height  of  5 m,  and  a wind  from  a 45°  azimuth  at 
10  m/sec.  Maximum  log^g  intensity  is  -2.781. 

11.  Intensity  isopleths  for  a clean  surface,  a source  elevation  of  30°, 

an  observer  height  of  150  km,  and  a wind  from  a 45°  azimuth  at 

10  m/sec.  Maximum  log^g  intensity  is  -4.661.  The  top  curve  is 
the  horizon. 

12.  Intensity  isopleths  for  a clean  surface,  a source  elevation  of  30°, 

an  observer  height  of  150  km,  and  a wind  from  a 90°  azimuth  at  10  m/sec. 

Maximum  log^g  intensity  is  -4.763.  The  top  curve  is  the  horizon. 

13.  The  same  as  Fig.  12,  except  that  an  oil  slick  is  added.  Maximum 
log^g  intensity  is  -4.366.  The  top  curve  is  the  horizon. 

14.  Intensity  isopleths  for  a clean  surface,  a source  elevation  of  30°, 

an  observer  height  of  150  km,  and  a wind  from  a 0°  azimuth  at  1 m/sec. 
Maximum  log^g  intensity  is  -4.112.  The  top  curve  is  the  horizon. 

15.  Intensity  isopleths  for  a clean  surface,  a source  elevation  of  30°, 

an  observer  height  of  150  km,  and  a wind  from  a 90°  azimuth  at  ^ 

1 m/sec.  Maximum  log^g  intensity  is  -4.101.  The  top  curve  is  the 
horizon. 

The  same  as  Fig.  14,  except  that  an  oil  slick  is  added.  Maximum 
log^Q  intensity  is  -3.985.  The  top  curve  is  the  horizon. 


Fig.  16. 
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Fig.  17.  The  same  as  Fig.  15,  except  that  an  oil  slick  is  added.  Maximum 
log^Q  intensity  is  -4.000.  The  top  curve  is  the  horizon. 

Fig.  18.  Intensity  isopleths  for  a source  displaced  100  m horizontally  from 
the  observer,  a source  height  of  62.74  m,  an  observer  height  of 
5 m,  a clean  surface,  and  a wind  from  a 45°  azimuth  at  10  m/sec. 
Maximum  log-jQ  intensity  is  -4.877. 

Fig.  19.  The  same  as  Fig.  18,  except  that  the  source  height  is  13.75  m. 
Maximum  log^g  intensity  is  -3.730. 

Fig.  20.  Isopleths  of  degree  of  polarization  (percent)  for  a source  displaced 
100  m horizontally  from  the  observer,  a source  height  of  62.74  m, 
an  observer  height  of  5 m,  and  a clean  surface. 

Fig.  21.  The  same  as  Fig.  20,  except  that  an  oil  slick  is  added. 
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